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Degenerate Fermi gas interacting with molecules near Feshbach resonance is unstable with respect 
to formation of a mixed state in which atoms and molecules coexist as a coherent superposition. 
Theory of this state is developed using a mapping to the Dicke model, treating molecular field in 
the single mode approximation. The results are accurate in the strong coupling regime relevant for 
current experimental efforts. The exact solution of the Dicke model is exploited to study stability, 
phase diagram, and nonadiabatic dynamics of molecular field in the mixed state. 



Feshbach resonance scattering [1-4], at which pairs of 
atoms can bind to form molecules at the same energy, 
has been used to demonstrate new coherence phenom- 
ena in cold atom systems. Those include, notably, the 
reversible coherent atom-molecule transitions [5,6] which 
can be accompanied by the Bose-Einstein condensation 
of molecules [7-9]. Recently, in search of fermionic con- 
densation, the focus shifted to Feshbach resonance in cold 
fermion systems [10-12]. 

The physics near the resonance in a macroscopic sys- 
tem is sensitive to the effects of quantum statistics. 
In particular, at positive detuning from the resonance 
molecules can coexist with fermions [13-16], stabilized 
by Pauli blocking of molecular decay into the states be- 
low the Fermi level. 

The stability and properties of the mixed state depend 
on the interaction effects. Below we argue that the in- 
teractions greatly enhance stability of the atom-molecule 
mixture, and lead to molecules and atom pairs hybridiz- 
ing to form a coherent state. We address the problem of 
molecules interacting with atoms by mapping it onto the 
Dicke problem [17] of two-level systems coupled to a Bose 
field. This problem, being exactly solvable [18], allows 
to describe the experimentally relevant regime of strong 
coupling. In the Feshbach resonance case, the two-level 
systems represent fermion pair states which can be occu- 
pied or empty, while the Bose field represents molecules. 

The coupling to molecules at positive detuning from 
Feshbach resonance enhances pairing interaction between 
fermions, which is expected to stimulate BCS superflu- 
idity [20,13-16,21,22]. In addition, as noted by Timmer- 
mans et al. [13] and Holland et al. [14], the strong cou- 
pling BCS condensation, with the critical temperature up 
to a fraction of Ep, may depend on the presence of molec- 
ular field. This conclusion was strengthened by a micro- 
scopic analysis carried out by Ohashi and Griffin [21], by 
Milstein et al. [22] , who refine the approach of Ref . [14] . 
Bruun and Pethick [16] studied noncondensed molecules 
coexisting with the Fermi gas at positive detuning using 
an effective theory of strong coupling formulated in terms 
of low energy parameters. It was noted that strong many- 
body effects exist even for detuning well above the Fermi 



energy. The important role of molecular field at positive 
detuning has been reemphasized recently by Falco and 
Stoof [15], who argued that a BCS-BEC crossover takes 
place in this region. 

In this article, we focus on the effects of molecule-atom 
hybridization and develop an approach allowing to han- 
dle this problem in the strong coupling regime. This is of 
interest, since the experiment deals with systems where 
the atom-molecule coupling, measured in the Fermi en- 
ergy units, is very large. We will see that molecule-atom 
mixing occurs in this situation in the range of detuning 
much larger than the Fermi energy, i.e. on the energy 
scale very different from that of fermionic condensate. 
The energy scale for the latter, set by the pairing inter- 
action strength, expected to reach 0.2 Ep at best [22], is 
much smaller than the atom-molecule interaction. Thus 
accurate results can be obtained with the help of a simple 
analysis which ignores direct pairing interaction between 
fermions, and relies on the exact solution of the atom- 
molecule dynamics. 

Below we analyze stability of fermions with respect to 
molecule formation, and obtain a phase digram. There 
is a fairly wide region around the resonance, span- 
ning both positive and negative detuning, were atoms 
and molecules coexist, forming a coherent state. At 
strong coupling, this region has width of the order of 
g^n/ Ep, a quantity which different estimates [15,16] put 
between few tens and few hundred Ep for current ex- 
periments [10-12]. Also, we exploit the Dicke problem 
solution to obtain nonlinear oscillations of the molecular 
field, in which population coherently oscillates between 
molecular and atomic components. The results of stabil- 
ity analysis are verified by comparing to the exact solu- 
tion and to the thermodynamic ground state properties. 

We consider the problem of a Fermi gas interacting 
with molecules in a single mode approximation which 
takes into account only the lowest energy molecular state: 
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with ap, Op 



+ 



.g ^ (6+Cp + h.c.) (1) 



and 6, h the atom and molecule opera- 
tors, a the fermion spin, and w the energy of a molecule. 
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The atom pair creation and annihilation operators Cp = 

^ (a-piapi + apia-pt)' 4 = 72 ("pT"-Pi + ''-pf4l) 
describe pairs of fermions in a spin singlet state that un- 
dergo conversion into molecules at Feshbach resonance. 
The approximation (1) is justified by the analysis below 
which finds that the energy gained by a formation of a 
mixed atom-molecule state, with all molecules occupying 
one state, is large compared to Ep, which allows to limit 
consideration to a single molecular state. 

The utility of the single mode approximation (1) is that 
it turns a difficult many-body problem into a well-known 
exactly solvable problem. The mapping is achieved by 
identifying the pair operators Cp, c+ with pseudospin 
Pauli operators [24] 

Cp = a-^l{a;-iay), = ^\{a; + ia^ (2) 

and noting that their product gives the particle number 
operator rip = a+Op in the subspacc of the many-body 
Hilbert space in which both states p and —p are simulta- 
neously filled or empty, 2c+Cp = rip + n-p = 0,2. More 
formally, defining Up = [(T+ , ] , one verifies that the 
standard Pauli spin commutation relations hold: 

[a+ , a;] = -2a+ , [a- , a|] = 2a; (3) 

This enables one to bring the Hamiltonian (1) to the form 
containing the spin variables only, 

^ = E' (^'^P + 9b^t + ab^^'i) + ^^^^ (4) 

where the sum is taken over singlet pair states with 
momenta p and —p. We note that the states with 
Up + ri-p = 1, with only one of the p and —p particle 
states filled and the other one empty, are decoupled and 
do not participate in the dynamics defined by (4). The 
reason for this decoupling is that these states have not 
enough particles to form a molecule, but also one particle 
too many to contribute to molecule dissociation. 

The spin-boson problem (4) is the Dicke model of 
quantum optics [17-19]. Hcpp and Lieb [18] found that 
the Hamiltonian (4) is integrable, and constructed exact 
many-body states. Besides the total particle number 

7V = 26+6 + ^(l + a^) (5) 

p 

there are also infinitely many nontrivial conserved quan- 
tities underpinning the exact solubility. 

The problem (4) resembles in many ways the BCS pair- 
ing problem. The latter is also integrable, which allows 
to obtain the full energy spectrum, and to construct non- 
trivial conserved quantities in a closed form [23]. In fact, 
the above pseudospin trick has its roots in the BCS prob- 
lem, where it was invented by Anderson [24] as an inter- 
pretation of Bogoliubov mean field theory. 



Here we employ the Hamiltonian (4) to assess stabil- 
ity of the Fermi gas with respect to molecule formation. 
The spin dynamics described by (4) is of the Bloch form, 
(T = i[H,(j] = 2hp X (7, with an effective magnetic field 
hp = {gb' ,gb" ,p'^/2m), where b = b' + ib" is a c-number 
describing the molecular state in macroscopic limit. 

The Bloch equations of motion for the spin components 

, (Tp and b take the form 

i(T+ = -(pV™)o-p + .9^0-p, i^p = ip'^/m)<Tp - gbap (6) 

i^p = '^gbup - 2gb*ap , ib = g^ a' + ub (7) 

p 

Prom a mathematical standpoint, Eqs.(6),(7) describe 

collective dynamics of a Bloch spin 1/2 ensemble, with 
the coupling between the spins provided by the 'magnetic 
field' hp transverse components which depend on the spin 
variables via an equation for b. Physically, the transverse 
spin components cjp characterize coherence between the 
filled and unfilled pair state, while describes the num- 
ber of pairs. 

Since the field 6 is a c-number, the operator equations 

(6), (7) are linear, and thus the spin operators expecta- 
tion values are subject to evolution equations of the form 
identical to (6), (7). In the absence of molecules, we have 
6 = 0, and all the spins are aligned in the ±z direc- 
tion, with probabilities determined by occupation of pair 
states: (crp = p-^ - p^ = - {1 - Up)^ = "^np - 1, 
where r?,p = {e^^'P /2'"-A') 1)"^ in thermal equilibrium. 
This state, containing only fermions but no molecules, 
(6) = (cr^) = 0, is stationary for the problem (6), (7). 

To assess stability with respect to molecule formation, 
we linearize Eqs.(6),(7), introducing 6ap,Sb oc e"*"^*. 
Sap ,5b* oc e'^**. Prom the coupled linear equations for 
Sa~ and 6b we obtain the eigenvalue equation 

^ = ^ + 9"Jl^P^ (8) 
p^/rn — A 

p ' 

To make the formally divergent sum over p well-behaved, 
following Ohashi and Griffin [25], we renormalize u) by 
subtracting the term 6ui = g^ ^p{p^ / 'm)~^ . The shift 
Lu —>■ Lu — Suj brings the position of the Feshbach reso- 
nance to w = for zero particle density, while Eq.(8) 
transforms to 

A = a;+g^vf l7~\ +^) (9) 

with the sum now converging at large p. 

The solution of Eq.(9) can be real or complex, depend- 
ing on the value of u). Complex-valued A = A' -|- iA" 
indicates an instability, with A" describing the instabil- 
ity growth rate. Numerical analysis of the solutions of 
Eq.(9) and simple analytic arguments reveal that (i) the 
real part A' is a monotonic function of io\ (ii) the instabil- 
ity occurs in an interval ujq < ui < lo\ with the threshold 
values Wo,! being a function of temperature. 
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The values wo.i can be inferred by noting that the com- 
plex A becomes real at w = wo.ij which gives the condi- 
tion A" — 0. When does Eq.(9) admit real solutions? 
This is possible for A < 0, as well as for A = 2fi, since 
2np — 1 changes sign at p = pp- (For all positive A ex- 
cept A = 2/i the residue (tTp) — 2np — 1 generates a finite 
imaginary part of A.) With A = 0, 2/i one obtains 

-o-VE4t^ (10) 
p 

2 f 1 ^ 2ri„ 1 \ 

wi = 2/i + .g > — ^-j— 

^—^ /m — p^ /m ) 

This indicates that atoms are stable at cj > wi, 
metastable at w < and at wq < ^ < t^i can ex- 
ist only in a state coherently mixed with the molecules 
(Fig. 1). We note that, since Wq < and Wi > 2/i, co- 
existence is favored by interaction. Moreover, at strong 
interaction, the detuning range where coexistence takes 
place becomes very large: Aw ~ g^n/Ep :s> Ep. 




Detuning co, |iK 



FIG. 1. Phase diagram of coupled atom-molecule system 
obtained from Eq.(lO) for ^'^K system [10] at particle density 
n ~ 1.8 X lO^^cm"^, Fermi energy Ep ~ 0.35/iK, and cou- 
pling strength g'^n/Ep ~ 60/iK. (The coupling was estimated 
using the microscopic theory of Feshbach resonance developed 
by Falco and Stoof [15], applied to the conditions of the JIIA 
experiment [10]). Inset: Effective potential schematic illus- 
trating the behavior in the three regions. 

The upper temperature at which cjq = is de- 
termined by the condition /i(T) = 0. For a two- 
species Fermi gas of total particle density n one has 
n = 2J2pnp{fi = 0) = 0.0972(m//3)3/2 which gives 
= 0.9885 i?_F. Interestingly, at low temperature 
T T^, the instability is pushed to higher detuning, 
uji — 2^ + (/^i^ln(/i/r), due to a BCS-hke log divergence 
at the Fermi level p — pp. 

It is instructive to look at the JILA experiment param- 
eters (Fig. 1). The estimate of coupling Auj — g^n/Ep ~ 



60/iK w 8 MHz gives a typical energy gained by the sys- 
tem via molecules and atom pairs hybridization, which 
is much larger than Ep. This leads to pair size in the 
mixed state ~ ?i/(2mAijj)^/^ much smaller than fermion 
wavelength Pp^. This indicates that the kinetic energy 
of atoms and molecules does not play a significant role, 
justifying the single mode approximation. 

Nonlinear dynamics at instability can be found with 
the help of the mapping to Bloch spins. Defining ~ 
{^{(Tp ± icTp)), Tp = (cTp), and rescahng gb 6, we write 

ir; ^ip^/m)r; ~br;, ir; ^ 2br; - 2b* (11) 

zb = ujb + g^J2'^P (12) 
p 

Since the norm is preserved by Bloch time evolution, 
Tpl^ = 4r~r+ -I- (r^)^ is conserved for each spin. We 
apply rotation, 

r;^e-^^%, r;^e"'\;, 6 ^ e'^^ (13) 

with the value rj to be determined later. This is equiva- 
lent to shifting p^ /ra ep ~ p'^ /m — r] and lu —> lu — i]. 

The resulting problem possesses real-valued solutions 
which can be obtained from the standard ansatz [26] 

r- = Apb + iBpb , = Dp - Cpb^ (14) 

Substituting this into Eqs.(ll),(12), from the real part 
of the equation or r~ and from the the equation or 
one finds the relations between the ansatz parameters 
Ap — epBp, Cp — 2Bp while the imaginary part of the 
equation or r~ generates a set of equations 

Bpb + epb - b{Dp - Cpb^) = (15) 

The constant of motion |rp|^ = r+ + (j'p)^ provides a 
first integral of Eq.(15): 

4 (ey + i?) + {2b' - Dp/Bp)^ = B;'\vpf (16) 

where we expressed Ap and Cp through Bp. 

Evidently, since the function b{t) is the same for all 
spins, the dependence on p has to drop out of Eq.(16), 
giving a single equation for 6 of the form 

b^ = {b' - b^_){bl - b^) , b-<b+ (17) 

which is possible with the following choice of constants: 

Dp/Bp - el ^b'_+bl, Dl-\Vp\^^AblblBl (18) 

These equations determine the modulus of Bp and Dp 
only. The sign has to be determined from initial condi- 
tions: sgn Bp = sgnDp = sgnr^. 

The solution of Eq. (17) is an elliptic function b{t) = 
dn(t, K^) with = 1 - b''_/b\ [27], oscillating period- 
ically between 6_ and b+. At 6_ <C 6+, the solution 
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is approximately given by a train of weakly overlapping 
solitons 



7 



cosh 7 — tn) ' 



tn. = rn 



(19) 



(Fig. 2), where each soliton in Eq.(19) is a solution of 
Eq. (17) with 6_ = 0,b+^ 7. 
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/ b /b =0.01 
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FIG. 2. Time evolution (14) of fermion pair amplitudes 
with different energies obtained for the molecular field of a 
soliton train form (19), shown in the inset. The Bloch sphere 
parameterization of pseudospin variables r^'^'^ = {o'p'^'^) is 
used. Pseudospins precess so that each state completes a full 
27r Rabi cycle per soliton. The red and blue curves correspond 
to the energies above and below the Fermi level. 

The quantities b± and the frequency rj must be deter- 
mined from the equation for b. One verifies that Eq.(12) 
is consistent with the ansatz (14) and obtains 



1 = 




Lu = r] 



Here = 2np — 1 corresponds to the energy distribution 
Up of fermions which depends on the initial state. For 
the parameters used in Fig. 1, by order of magnitude we 
estimate 7,t~^ ~ g^n/Ep w 8 MHz. This is much faster 
than typical energy relaxation rates, which justifies ig- 
noring relaxation effects in the dynamics. 

The properties at equilibrium can be understood by 
considering the limit b^ ^ b^ — bo, when oscillations are 
absent. The energy distribution Up can be easily obtained 
in the pseudospin picture, taking into account that each 
spin is presented with a tilted field hp = (60, 0,p^/2m— /x) 



which gives Up = 1/ (l -I- e'^'*'''!). The molecular field 60 
in the ground state is determined by 



E 



sgnep (1 - 2np 



Abl 



(21) 



along with the constraint N = 2&g/ + 2 J^p "•p- 

Here we use Eq.(21) to verify the above stability anal- 
ysis. To determine when the atoms can be stable with 
respect to hybridizing with molecules, we set 60 = and 
immediately recover Eq.(9) for the instability exponent 
A. The difference, however, is that 77 defined by Eq.(21) 
is real, while A is complex. Atoms' stability is thus in- 
deed equivalent to the existence of a real- valued solution 
of Eq.(9). One possibility to have such a solution is to 
set 77 — 2/i, which eliminates the log divergence in (21) 
at Ep — 0. The other possibility is to have /i, 77 < 0. Put 
together with the properties of equilibrium state at finite 
60, this confirms the above estimate of the coexistence 
region (10) and the conclusion that pure atom state is 
metastable at the detuning lu < ujq. 

In summary, this work provides a phase digram and 
an exact solution for the atom-molecule dynamics in the 
regime of strong coupling. The characteristic energy 
scales are estimated to be much larger than Ep, which 
makes the Dicke model approximation ignoring molecular 
dispersion as well as the BCS fermion pairing effects, ac- 
curate enough. A wide atom-molecule coexistence region 
is predicted in which atom pairs and molecules hybridize 
into objects of size much less than Fermi wavelength Pp^ . 

After having completed this work we became aware 
of the article by Andreev, Gurarie and Radzihovsky [28] 
which exploits the mapping to the Dicke model, while 
focusing on the weak coupling limit. 
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